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Abstract 

We study the force vs extension behaviour of a helical spring made of a thin 
torsionally-stiff anisotropic elastic rod. Our focus is on springs of very low heli- 
cal pitch. For certain parameters of the problem such a spring is found not to 
unwind when pulled but rather to form hockles that pop-out one by one and lead to 
a highly non-monotonic force-extension curve. Between abrupt loop pop-outs this 
curve is well described by the planar elastica whose relevant solutions are classi- 
fied. Our results may be relevant for tightly coiled nanosprings in future micro- and 
nano(electro)mechanical devices. 

Key words: thin anisotropic elastic rod, helical spring, cascade unlooping, planar 
elastica, force-extension response, loop pop-out, hysteresis 



1 Introduction 



Take an elastic rod or strip that is coiled when left free and pull the ends apart 
without fixing their orientation. We can do this by attaching strings to the 
ends and applying a tensile force to the strings (see Fig. [T]). If we imagine the 
rod to be held between two parallel glass plates, forcing it to adopt an essen- 
tially planar configuration, then we would observe the formation of a number 
of loops given by the number of initial coils. The loops have nowhere to go 
under increased stretching of the rod and grow into tight hockles with a highly 
localised curvature. Without the constraining walls the planar configuration 
becomes unstable at a certain force and unlooping occurs. However, this un- 
looping force can be very high, especially if the rod has a flat cross-section, 
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Fig. 1. Photographs of tensile tests on a plastic helical strip of small pitch. The 
unstressed strip has four loops, which are successively lost as the force is increased. 



and the structure may get damaged (e.g., through plastic deformation) before 
loop pop-out occurs. Indeed, hockles (or kinks or snarls) are a serious concern 
not only in the use of garden hoses but also in such slender industrial struc- 
tures a s marine pipe lines, mooring ropes, textile yarns and communication 
cables flCovnel . Il990h . Pipes, ropes, yarns and cables often acquire intrinsic 
curvature due to creep from storage on a drum, putting them at risk for the 
above hockling scenario. 



Helical springs have also recently been fabricated at t he micro and nano scale 
flGao et all 120051 : IChol . 120061 : IPrinz and Golodl . l2006f ). where they hold great 
promise as essential parts of future micro- or nano (electro) mechanical devices 
such as che mical or bi o logica l sensors, cantilevers, resonators, inductors, actu- 
ators, etc. (jBell et all 120061 ). Nanohelices have been synthesised from differ- 
ent materials leading to structures with unusual mechanical properties, and 
force-extension experime nts on nano s prings , nanocoils and n a nobel ts are ac- 
tively being performed ( Chen et all 120031 ; iKratochvil et all 120071 ). For in- 
stance, ZnO and Si 3 N 4 springs have been shown to be superelastic, i.e., they 
can be stretched to an a lmost straight line and when releas ed they restore their 
shape without damage (jGao et all 120061 ; ICao et all 120081 ). An i nteresting de- 



velop ment is the fabrication of nanosprings with very small pitch (jZhang et al. 



20061 ). Such springs are particularly desirable because they allow for a large 
magnetic flux density. 
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Fig. 2. Schematic diagram of a pulled spring of helical angle 8. A flat strip is drawn 
to emphasise the anisotropic bending stiffness of the rod. 

In this paper we consider the problem of successive hockling and unloop- 
ing of a helically coiled stri p of low pitch under tensile load. Th is problem 
was briefly touched upon in (jStarostin and van der Heijdenl . 120081 ). where we 
modelled the strip as an inextensi ble helical shell mo tivated by multistabil- 
ity observed in cholesterol ribbons (ISmith et all l200ll ) . Here we do not make 



the inextensibility assumption an d model the strip i n the more common way 



as a transversely anisotropic rod (jZhou et all 120051 ). The undeformed rod is 



curved about the axis of least bending stiffness (so that our spring is close to a 
clock spring, rather than, say, a slinky toy; see Fig. [2]). It is worth pointing out 
here that although the rod will be assumed to have an inextensbile centreline, 
this rod model is not equivalent t o the inextensible strip model i n the limit 
of infinitesimal width, as shown by lStarostin and van der Heijdenl (120071 ). The 
two models apply to different circumstances. 



We choose zero-moment boundary conditions because these seem to be a rea- 
sonable choice to avoid kinking as the rod is free to unwind (the loops are not 
'locked in' as in the case of clamped boundary conditions). Nevertheless, we 
show that under certain conditions on the elastic and geometrical parameters 
a tightly coiled strip when pulled behaves in a highly nonlinear way. Under 
small forces the spring's force-extension behaviour is Hookean, but at larger 
extensions the loading curve develops hysteresis cycles as more and more loops 
pop out. These pop-outs form brief and abrupt three-dimensional excursions 
in a process that is otherwise remarkably well approximated by planar elastic 
behaviour. 



This planar behaviour is described by the Euler elastica subject to a constant 
intrinsic curvature. The fact that a low-pitch helical spring is found to jump 
between planar states justifies a revisit of the Euler elastica and we give what 
appears to be a new classification of its solutions subject to zero-moment end 
conditions. The corresponding force-extension curves provide the scaffold of 
the tensile response of the helical spring. A t larger pitch no tight hoc kles occur, 
although multistability may be observed (IKessler and Rabinl . 120031 ) . 
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Multi- l oop ho ckling and unlooping does not se em to have be e n stud ied before. 
Coyne only considered a single loop. lYabuta et al.l (119821 ) discussed 



cable loop stability assuming a single tightly coiled helical turn. Both works 
assume the unstressed rod to be straight and avoid boundary conditions by 
assuming the cable to be infinitely long. We model the strip as a finite-length 
intrinsically curved elastic rod of anisotropic (flat) cross-section. 

The organisation of the paper is as follows. In Section 2 we set up the equi- 
librium equations for elastic rod statics. In Section 3 we classify the planar 
solutions of this strip and find the force-extension curves which the curves for 
a coiled strip will be compared against (in Section 4). The main result of the 
paper is Fig. [TDl which shows the hysteretic force-extension behaviour char- 
acterising multi-loop pop-out. We briefly contrast this behaviour with that of 
springs of different geometric and elastic parameters. Section 5 concludes the 
paper. 



2 Anisotropic rod model 

We consider a thin inextensible elastic rod with non-circular cross-section and 
free from distributed external loads. Let the centreline of the rod be r(s), 
where s G [0, 1] is arclength and the length of the rod has been scaled to 1. 

The two principal axes of the cross-section and the tangent to the centreline, 
t — r' , constitute an orthonormal material frame {d\, d 2 , d^}, where d% = t 
and a prime denotes differentiation with respect to s. By orthonormality we 
have the kinematic equations 

d'^uxdi, i = l,2,3, (1) 

where u = {u\, u 2 , us) T is the strain vector written in the material frame. 

Let F(s) and M(s) be the resultant internal force and moment that the 
material of [s, si) in the rod exerts on the material of (s , s], where s < s < S\. 
In case there a re no distrib uted loads the force and moment balance equations 



take the form (JLoveJ, 119271 ) 



F' — 0, M' + t x F = 0. (2) 

We close the system of equations by specifying (linear) constitutive relations 
Md 1 = A(ui-mio), 

Md 2 = B(u 2 -u 20 ), (3) 
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M -d 3 = C(u 3 -u 30 ), 



where we have allowed for non-zero curvature and torsion of the rod in its 
unstressed state. Here A and B are the principal bending stiffnesses and C is 
the torsional stiffness, while U\q and U20 are the intrinsic curvatures and M30 
is the intrinsic torsion. 

To these equations we apply the following zero-moment boundary conditions: 



x(0) 


= 0, 


x(l) = 


0. 


y(o) 


= 0, 


y(i) = 


0. 


s(0) 


= 0, 


^(i) = 


F z , 


M x (0) 


= 0, 


M X {1) -- 


= 0, 


M y ((0) 


= 0, 


My(l) - 


= 0, 


M z (0) 


= 0, 


di(l)-* 


= 



(4) 



where we have introduced components as follows: r = xi + yj + zk, F = 
F x i+F y j+F z k, M = M x i+M y j+M z k, with {i, j, k} a fixed laboratory frame 
that has its origin at one end of the rod and k aligned with the (constant) force 
in the rod. Thus we impose the constraint that both ends of the rod remain 
on the line through k. Together with the six orthonormality conditions 



d»(0) • dj(0) 



j> z = l,2,3, 



(5) 



where 6ij is the usual Kronecker symbol, Eq. (TJJ gives a set of 18 boundary 
conditions for the 18 equations for the components of r, F, M and dj (i = 
1,2,3). 



We solve this boun dary- value problem using the continuation code AUTO 
(IDoedel et all I1998T 1 , which allows for the tracking of solutions as parameters 
are varied and also detects bifurcations. 



If we strongly penalise bending about the first material axis by pushing the 
ratio A IB to infinity then the m aterial frame will be locked onto the Frenet 
frame (IKessler and Rabinl. 120031). The corresponding equilibrium equations 
were derived by IStarostin and van der Heijdenl (120071 ) . 



Before solving the equations, in the next section we consider the special case 
where the unstressed rod is a multi-covered ring and the deformations are 
planar, described by Euler's elastica. 
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3 Planar elastica with intrinsic curvature 



3.1 Solving the equilibrium equations 



In this section we consider the planar shapes of a twist-free thin inextensible 
elastic rod. This degenerate planar case is important for the following anal- 
ysis because it will help us better understand the behaviour of the low-pitch 
spring under tension. When left free of external forces and moments, the rod 
is assumed to take on the shape of a multicovered ring of radius R . Let it 
wind n times so that the intrinsic curvature becomes xq = 1/Rq = Inn (recall 
that the rod has unit length). Let s\ and S2 be the arclength coordinates of 
the ends of the rod so that s 2 = Si + 1. 

We are interested in the equilibrium configurations under an applied uniaxial 
non- vanishing tensile force F at the ends of the rod. Equations (j2j) in this case 
reduce to one equation for the tangent force F t = F ■ t: 

{F t 'f = {F t + H){F 2 -F 2 \ (6) 



where if is a constant of integration and we have set B = 1. Recall that the 
force vector is constant, hence F 2 = const. The intrinsic curvatur e x n does no t 
enter the equation, which is therefore that of the Euler elastica (iLovd . 119271 ). 



By introducing the angle d between the fixed direction of the force and the 
tangent to the rod we can write F t = F7 = F cos d and Eq. ([6]) takes the form 



( 7 r = (i/ + F 7 )(i-y) 



(7) 



(cf. Eq. (2.8) in (jStarostinl . 12004 )). Integration of this equation gives 

7 = l-2sn 2 (fts), (8) 



where Q 2 = (H + F)/4. For non-inflexional solutions the elliptic modulus k is 
given by k 2 = 2F/(H + F), where we require that H > F. For the force we 
can thus write F = 2k 2 Q 2 . 

We introduce planar coordinates (x, z) with z chosen along the direction of 
the force. We may then write x' = ±\/l — 7 2 , z' = 7, with 7 given by Eq. (jSJ), 
and the centreline is found by integrating these equations (or they can be 
immediately obtained from Ilyukhin's equations reduced to the planar case 
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( iStarostin 



x = 



(9) 



where E(w, fc) = / " dn 2 w du> is the incomplete elliptic integral of the second 
kind. 



The curvature of the centreline is x = k 2 Q 2 x (we note that always x > 0). 

The shape described by Eq. (|9]) is periodic with period 2K(k)/Q, where K(k) 
is the complete elliptic integral of the first kind. Moreover, the centreline is 
symmetric with respect to the x-axis so that x(s) = x(—s) and z(s) = —z(—s). 
The curvature monotonically decreases on the intervals (2mK(k)/Q, (2m + 
l)K(k)/Q), m G Z, and monotonically increases on all other intervals. 

We apply moment-free boundary conditions, i.e., = x for i = 1, 2, which 
gives 



This implies that x(si) = x(s2) and that the force acts along the end-to-end 
line. 



3.2 Classification of solutions 
3.2.1 Primary symmetric solutions 

We first consider symmetric shapes such that the midpoint of the centreline 
coincides with the extremum of the curvature, i.e., (si + S2}/2 = mK(k)/fl, 
m 6 Z; it is a maximum for even m and a minimum for odd m. Then the 
condition s 2 — Si — 1 of fixed length transforms into Sj = mK(k)/Q — (— 1)*/2, 
i = 1,2. Substitution of the latter into Eq. (flOl) results in one of the two 
following conditions, depending on the parity of m: 



Tin = Q dn(Qsi) 



i = 1,2. 



(10) 



Tin = Qdn (fi/2) for m even, 



Tin = VlVl — k 2 nd (Vl/2) for m odd. 



(11) 
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Fig. 3. Elastica configurations for n = 1. Left: f-shape solution. Right: a-shape 
solution. The end curvature in both cases is = 2nn. The elliptic moduli of the 
two solutions are the same. The force acts horizontally. 

The end-to-end distance Az = z(si) — z(s2) is computed from the second 
expression in Eq. ([9]) as 



For given k (and m) we consider Eq. (fTTj) as an equation for Q. For k = 0, 
each of Eq. fill I) has the unique solution Q = im. By continuity in k we can 
conclude that for given n there exist only two branches of symmetric solutions 
through the undeformed initial state (one for even m and one for odd m). 
Depending on the parity of n — m these two types of solutions have different 
shapes: if n — m is even, then we call them v -solutions, otherwise we call 
them ot-solutions (see Fig. [3]). (This labelling is prompted by the superficial 
resemblance of the centreline configurations to the Greek letters.) These two 
types may be thought of as opened (v) and tightened (a) coils. 

Equation ffTTj) is solved numerically for increasing values of k. Force-extension 
curves for both types of symmetric solutions are shown in Fig. HI Note that in 
this and all the subsequent figures we always show the negative external force 
that balances the positive internal force in the rod. The force is furthermore 
scaled by (27m) 2 with n = 4 for easy comparison with numerical results for 
the 4-turn helical spring in Section 4. 

For the v-shapes we can approximate the solution for small fc < 1 as O = 
irn + y^tt 3 ™ 3 *; 6 + 0(k s ). Then F = 2k 2 Vl 2 = 2ix 2 n 2 k 2 + 0(k 8 ), while from 
Eq. (|T2|) for the extension Az we can write 



For the a-shapes the approximate solution of Eq. ( TTTT) is different: Q = 




(12) 




(13) 




+ 
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Fig. 4. Force-extension curves for the rod with intrinsic curvature xo = 2im. Sym- 
metric branches passing through the origin are shown for n = 1, ... ,8 (increasing 
from right to left). Solid curves are for u-solutions, dashed curves for a-solutions. 



k A ) + 0{k 8 ) and 



457 



Az = -k 2 + —k A + 

8 64 \ 8192 



7r 2 n 2 ~ 



96 



k 6 + 0{k 8 ). 



(14) 



These approximations show that for small force and extension the slope of 
the force-extension curve of both types of solution is the same, i.e., the linear 
response is F = y7r 2 ?2 2 A2;. Indeed the curve for one type of solution can 
be considered as a continuation of the curve of the other type of solution in 
the other direction (negative Az and positive F z ). For larger extensions the 
two curves diverge. For relatively large forces the curve for the a-solution 
corresponding to n approaches the curve for the ^-solution corresponding to 
n + 1, as illustrated in Fig. HI Other symmetric solutions emerge at finite 
extension and force (resp., k), as we shall see below. 



3.2.2 N on- symmetric and secondary symmetric solutions 

There are also non-symmetric configurations (we call them d-shapes) such 
that the length of the rod comprises an integer number of periods j, i.e., 
2K(k)j/Q = S2 — Si = 1, but the midpoint is not at maximum or minimum 
curvature (see Fig. Hence, Q = 2K(k)j and the coordinate for the initial 
point of the rod Si can be found from the condition x(si) = x = 27m, which 
gives 

nn = 2K{k)j dn{2K{k)js 1 ). (15) 
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Fig. 5. Example of a non-symmetric $-shape solution for n = 3 and j = 3. Both 
end curvatures are >cq = 2nn. The force acts horizontally. 

The extension of non-symmetric solutions is computed from Eq. (Q with Q = 
2K(k)j: 



where E(/c) is the complete elliptic integral of the second kind. The force is 
given by F = 8k 2 K(k) 2 j 2 . Force vs. extension curves are shown in Fig. [6] for 
./ 1 7. 

Note that the loading curves in Fig.[6]do not themselves depend on the intrinsic 
curvature, i.e., on the number of initial coils n, but their parametrisation by 
k does, as is seen in Eq. (1151) . For k = this equation has a solution only for 
j = n (then Q = irn) and, by continuity, this is also true for small k. For such 
k we can write, from Eq. ( fl~6l) . 

Az = k 2 (- + —k 2 + —k A + 0(k 6 )) . (17) 
V8 16 1024 v V y J 
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Developing the force into the series in k, we have, with j = n, 

F = n 2 n 2 k 2 (2 + k 2 + l -h 4 + 0(fc 6 )) . (18) 

Combining both expressions we find the linear approximation F = 16ir 2 n 2 Az, 
valid for small deformations. We see that the non-symmetric solution requires 
a force three times that of the symmetric one. 

As k increases, other roots of Eq. ( |T5l) appear. A new root comes into existence 
each time the Jacobi elliptic function dn reaches a maximum or minimum, i.e., 
for s\ = ;p (even m corresponds to a maximum, odd m to a minimum). A 

maximum can only occur for j < n, because j = 2K j fc - ) > 1 for k > 0, while a 
minimum can only occur for j > n, because ™ = 2K( - fc - ) ^ 1 ~ — < 1 for k > 0. 

j <j i j iff 

Thus, for given actual number of periods j the corresponding curve in Fig. [6] 
represents the shapes for various intrinsic curvatures. For n — j the entire 
curve is realisable; otherwise, only part of it is realisable with some offset 
from the origin. The offset increases with the difference between the number 
of initial coils n and actual periods j. 

The values of k where new roots emerge correspond to configurations at the 
intersection of the non-symmetric branch with secondary symmetric branches 
that do not pass through the origin (see Fig. [7]). For j < n, they are ^-shapes 
with their ends at points of maximum curvature, while for j > n they are 
a-shapes with their ends at points of minimum curvature (see Fig. [S]). In both 
cases the end tangent vectors are along the end-to-end distance and the z-axis. 

If we fix the intrinsic curvature, i.e., the number n, then the set of the non- 
symmetric solutions consists of 1) one complete curve j = n that begins in 
the origin and 2) parts of all other curves starting from the bifurcation points. 
There is only one bifurcation point on each j-curve and the symmetric branch 
is born here. The solutions on these branches are computed from Eqs. (TTTT) 
and (jT2l) with m = j for j < n and m = j + 1 for j > n. Fig. [9] shows 
non-symmetric and bifurcating symmetric branches for n = 4. It is conve- 
nient to label the symmetric branches with the number j. As the applied force 
increases, the symmetric j-branches for j > n (a-shapes) converge to non- 
symmetric branches with the same j, while the branches for j < n (-u-shapes) 
approach their two neighbouring non-symmetric curves with j±l, i.e., they ei- 
ther lose or gain one coil. Meanwhile, the non-symmetric j-branch approaches 
the primary symmetric n-branch with n — j + 1. 

Physically, solutions on the bifurcating symmetric branches, which do not go 
through the origin, can only be obtained by uncoiling the n-loop and recoiling 
into a shape with j loops. They cannot be held at zero force. 
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Fig. 7. Force-extension curves for a rod with intrinsic curvature xq = 2irn with 
n = 1, . . . , 10. Non-symmetric and some symmetric branches that do not pass 
through the origin are shown for j = 1,...,7. Bifurcation points are shown for 
various values of n. 

Summing up, we see that for given intrinsic curvature there exist 1) only 
three branches (symmetric a- and ■u-shapes and a non-symmetric ^-shape) 
that contain the relaxed state, 2) n — 1 associated pairs of opened and 
^-solutions, 3) a countably infinite number of tightened and a-soiutions. 
Clearly, no branch of opened or tightened shapes can approach the relaxed 
state of the n-covered ring. 



4 Pulling a helical spring — Numerical results 

Having classified all planar solutions we now turn to the three-dimensional 
deformations of our helical spring. 

(ill) 

Fig. 8. Elastica configurations at intersection points of symmetric and non-symmet- 
ric branches for n = 4. Left: v-sh&pe solution for j = 3. Right: a-shape solution for 
j = 5. The force acts horizontally. 
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Figure [H)] gives the force-extension curve for an unstressed helix of unit length 
with 4 turns (i.e., n = 4) and helical angle 9 = 89.8° (i.e., a pitch an- 
gle of 0.2°), subject to the above boundary conditions. This helix has ra- 
dius i?o — sin #o/ (27rn), curvature x = 27m sin # and geometrical torsion 
t = 27m cos 8q. The elastic constants are fixed such that A/B = 10 and 
C/B = 1.7316, and we set u w = 0, w 2 o — x o an d u 30 = r . These values imply 
that the smaller bending stiffness is along the axis of intrinsic curvature. The 
value of C/B corresponds to a Poisson ratio of v — 0.05 if we assume the 
elastic constants to be at the boundary of the admissible region for homoge- 
neous cro s s-sect ions of arbitrary geometry, as given by the Nikolai inequality 



flllvukhinl . Il979f ) 



^ (ihJkIhb) - (19) 



(Here equality holds only for elliptical cross-sections.) All forces are scaled by 
{2im) 2 B. 

The figure shows that after a short stretch, during which the helix turns and 
shears, the force-extension curve very closely follows the (primary symmetric) 
planar elastica curve until the elastica's out-of-plane instability is approached. 
The curve then sharply veers away from the j = 4 elastica curve and solutions 
become three-dimensional. This is illustrated by the shapes depicted in Fig.[TTJ 
which bear good resemblance to the (identically labelled) actual strip shapes 
in Fig. [TJ After this excursion the curve lands on the (secondary symmetric) 
j — 3 branch, its solutions having shed one loop. Two more such unlooping 
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0.2 0.4 0.6 0.8 1 
Az 

Fig. 10. Force-extension curve for a low-pitch helix with n = 4 (thick curve) and 
approximating planar elastica curves. The diamonds indicate points of out-of-plane 
instability of the elastica. The enlargement reveals a sharp transition from 3D to 
2D behaviour. Labels 'a' to 'h' refer to solutions displayed in Fig. [11] (# = 89.8°, 
A/B = 10, v = 0.05.) 




h 



Fig. 11. Solution shapes along the curve in Fig. [TU1 Labels loosely correspond 
to those in Fig. [TJ Note the shearing behaviour at small force (inset), (n = 4, 
6 = 89.8°.) 

manoeuvres occur while the force-extension curve jumps from the j — 3 to 
the j = 2 and finally to the j '• = 1 branch. Associated with these jumps is 
hysteresis behaviour, i.e, the jumps occur at different points if the extension 
Az is decreased rather than increased. It is worth noting that all elastica curves 
visited by the spring curve are u-branches. The small pitch angle acts as an 
imperfection causing a rounding off of the elastica bifurcation diagram which 
would have a branch of three-dimensional solutions connecting the pitchfork 
bifurcations indicated by diamonds. Note that the location of these diamonds, 
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Fig. 12. Force-extension curves for different helical angles: 9q 
89° and 88° (n = 4, A/B = 10, v = 0.05). 



89.8° (solid), 89.5°, 



Ph 




Fig. 13. Force-extension curves for different anisotropics: A/B = 10 (solid), 8, 6 and 
4. Also shown are curves for isotropic (A/B = 1, i) and slinky (A/B = 1/10, ii) 
rods (solid), together with approximations (dashed) obtained by assuming an exact 
helical shape throughout, (n = 4, 9q = 89.8°, v = 0.05.) 

but not the curves, depends on the torsional stiffness C. Clearly, for the nearly 
planar solutions not to have self-contact, which is not modelled, we have to 
assume the strip to be sufficiently narrow. 



Figures [12] and [13] show the dependence of this unlooping behaviour on 9 and 
anisotropy A/B. It is seen that the unlooping cascade requires, in addition to 
both a small pitch and a small u, a relatively large anisotropy, otherwise the 
force-extension behaviour resembles that of a helix, which gradually unwinds 
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Fig. 14. Force-extension curves (solid) for a helical spring at 8q = 88° 
(Azo = cos #o = 0.03490) and 60° (Azq = cos#o = 0.5) compared against curves 
(dashed) obtained by assuming an exact helical shape throughout (n = 4, A/B = 10, 
v = 0.05). 

under an applied tension. To illustrate this, in Fig. [14] we compare force- 
extension curves for 6q — 88° and 60° with curves obtained by assuming an 
exact helical shape for the rod throughout the tension test. For this case, which 
would require special boundary conditions, exact expressions can be derived. 
These are given in the Appendix. We note that for a pitch angle as small as 
2° (#o — 88°) the agreement with the helical curve is already reasonably good, 
although the latter fails to account for multiple local minima. At relatively 
large pitch angles (e.g., #0 = 60°) the hysteresis behaviour disappears. 

For comparison, in Fig.[13]we also present force-extension curves for the slinky 
spring, whose axis of least bending stiffness is at right angles to the axis of 
intrinsic curvature (A/B = 1/10), as well as for the isotropic rod (A/B = 1). 
In both cases no hysteresis is found and the curves are well approximated by 
helical curves (included in dashed lines) . Slinky shapes are shown in Fig. [15] 
revealing that the slinky spring is pulled out into an almost perfectly heli- 
cal shape. Interestingly, rather than unwinding, the spring overwinds under 
increasing force. 



5 Concluding remarks 

We have analysed the problem of successive unlooping of a helical spring. 
Springs of sufficiently low pitch made of material of sufficiently large ratio of 
torsional to bending stiffness are found not simply to unwind when pulled, 
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Fig. 15. Slinky solution shapes along the curve labelled ii in Fig.[T3J The inset shows 
shape 'a' for a narrower strip to reveal the characteristic slinky response, (n = 4, 
O = 89.8°.) 



but rather to form hockles that successively pop out. The force-extension 
curves are non-monotonic with large parts well approximated by curves for 
the planar Euler elastica and connected by brief three-dimensional excursions. 
Similar complicated force-extension behaviour was found in the helical shell 
model proposed by IStarostin and van der Heijdenl (120081 ). In that study we 
showed that the spring force-extension curves can be approximated by (pri- 
mary) elastica curves (analogous to those in Fig. H]). In the present paper we 
have shown that the force-extension curves are even better approximated by 
the secondary elastica curves that bifurcate from non-symmetric branches (as 
in Fig. ED. 

A low-pitch helix can be interpreted as a 'perturbation' of a planar circu- 
lar spring. Our analysis gives an unfolding of the bifurcation diagram of the 
perfect planar elastica. The unlooping scenario in the unfolding describes the 
tensile response of an imperfect planar spring. 

Our study has been restricted to (quasi-) statical solutions. In practice the 
three-dimensional excurs ions will be rapid a nd would require dynamical mod- 
elling, as for instance in ( Goyal et all 2005 ). Our work also paid no attention 
to self-contacts. This is justified by the fact that no shapes we computed cross 
themselves (the only exception is the slinky shape at very small force). Thus 
choosing cross-sections small enough excludes self-contact. 



By choosing our torsional stiffness to satisfy the Nikolai condition as an equal- 
ity we take the maximum possible value for any solid cross-section, given the 
bending stiffnesses A and B and the Poisson ratio v. Even so, we need quite 
a small value of v to observe the unlooping cascade described in the previ- 
ous sections. However, moving away from solid cross-sections, materials (e.g., 
biopolymers) are now well-known that have a relatively large torsional stiffness 
corresponding to a small or even negative effective Pois son ratio. For insta nce, 
double-stranded DNA is exceptionally stiff in torsion ( Oroszi et all 2006 ). It 
is also worth pointing out that the inextensible helical shell model considered 
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by IStarostin and van der Heijdenl (120081 ) would exhibit the unlooping cascade 
for a larger range of v values because of its higher intrinsic torsional stiffness. 



Our results may be useful for the design of mechanisms for the deployment 
of ribbon-like structures. The nonlinear properties of the elastic response of 
helical strips could be exploited in mechanical sensors and actuators at the 
micro- or nanoscale, where low-pitch springs are being considered. 



A Force-extension curves for helical solutions 



Helical curves are characterised by constant curvature and torsion. Since the 
stiffness ratio A/B is generally taken relatively large in this paper (A/B = 10 
in Figs. [TU] and [T41) bending of the rod will predominantly be about d 2 , the 
direction of intrinsic curvature. Thus we seek helical solutions of the equations 
in Section 2 with Ui — 0, u 2 = >c = const., Us — t — const. It is convenient 
to introduce the helical angle 9 (cf. Fig. [2]) such that tan# = x/r. Force and 
extension for a unit-length helix are then parametrised in terms of 9 as follows: 



BC sin 2 6*0(5 sin 9 sin 9 + C cos #0 sin 9) sin(0 - 9 ) 
~ Rl sm9{Bsm 2 9 + Ccos 2 9) 2 ' ( 

Az = cos9, (A.2) 

where Rq = sin 2 o / x o is the radius of the undeformed helix and z is a co- 
ordinate along its axis. Note that the force is now pointing along the axis of 
the helix and not necessarily along the end-to-end vector. The radius of the 
deformed helix is given by 

R = Rq *™ 9(B sm 2 9 + Ccos 2 9) 

sin 9q (B sin 9 sin 9q + C cos 9 cos 9q ) 



In Fig. [TU force-extension curves for these helical solutions are included for 
comparison with the force-extension curves for a pulled helix. 
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